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Abstract 

We consider the spectrum of BPS saturated states in = 2 gauge theories in four 
dimensions. This spectrum may be discontinuous across real codimension one submanifolds 
of marginal stability in the moduli space of vacua. An example, which can be treated with 
semiclassical methods in the weak coupling limit, is the decay of quark-soliton bound states. 
For a quark and a soliton of electric-magnetic charge vectors Q and Q' respectively, we find 
that as the manifold of marginal stability is crossed, the number of soliton states changes 
by a factor of 2^'^ , where the dot denotes the symplectic product. 
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1. Introduction 

An important feature of = 2 supersymmetric theories in four space-time dimensions 
is the existence of a central charge Z in the supersymmetry algebra |jl|]: 

{Qa,Qf} = ea/.e^''^- (1-1) 

The central charge is in general a linear combination of conserved abelian charges with 
coefficients that depend on the vacuum moduli and the parameters of the Lagrangian. In 
a unitary representation of supersymmetry, the mass m of a state is bounded by 

m>\Z\. (1.2) 

If we limit ourselves to maximum spin 1 (i.e. we do not consider supergravity theories) 



we have the following representations of the supersymmetry algebra ([T^): When the 
inequality (|1.2| ) is not saturated, an irreducible representation contains five scalar, four 
spinor and one vector particle, i.e. a total of sixteen helicity states. For a so called BPS 
representation, which saturates ( |1.2| ), there are two possibilities; either two scalar and one 
spinor particle, or one scalar, two spinor and one vector particle, for a total of four or eight 
helicity states respectively. 

Particles in a BPS representation enjoy the following stability property: For a BPS 
state with central charge Z and mass m to decay into states with central charges Zx and Z2 
and masses mi and m2, we must have Z = Zi + Z2 by charge conservation. But this means 
that m = \Z\ < \Zi\ + \Z2\ < mi + m2 by the triangle inequality, so this decay cannot take 
place. These inequalities are saturated exactly when the would-be decay products are also 
BPS states and the central charges Zi and Z2 are aligned in the complex plane. This will 
happen on a real codimension one subvariety of 'marginal stability' in the moduli space 
of vacua. As explained in 0, the spectrum of BPS one-particle states must sometimes be 
discontinuous across such a boundary. This phenomenon is well understood in two space- 
time dimensions where the number of particles appearing or disappearing in this way 
can be calculated from Picard-Lefschetz theory in the case of Landau-Ginzburg theories, 
or from the Tt{—l)^Fexp{—f3H) index and the tt* topological-anti-topological fusion 
equations for general N = 2 theories. 

The purpose of this article is to study such discontinuities in the context of A = 2 
gauge theories in four dimensions. In the simplest example, the A = 2 SU{2) gauge 
theory without extra matter there are only two conserved abelian charges, namely the 
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electric and magnetic charge for the unbroken ^7(1) gauge group. Consequently, there is a 
single curve of marginal stability where their coefficients in the central charge are aligned. 
Unfortunately, this curve lies entirely in the strong coupling regime , where it is difficult 
to exhibit the discontinuity of the spectrum explicitly. In more complicated examples, i.e. 
with larger gauge groups and/or extra matter, there are more conserved abelian charges 
and the manifold of marginal stability depends on the quantum numbers of the particles 
involved in the decay. (The set of all such curves, regardless of which particles actually 
exist in the spectrum, is in fact dense in the moduli space.) More important for our 
purposes is that models with extra matter have dimensionful parameters, namely the bare 
masses of the matter fields. When such a parameter is large compared to the dynamically 
generated scale of the theory, marginal stability might occur also at weak coupling, where 
semiclassical methods are applicable. 

This paper is organized as follows: In section two, we quickly review some aspects of 
N = 2 gauge theories, and in section three, we discuss their BPS spectra at weak coupling. 
In section four, we calculate the number of quark-soliton bound states that disappear as 
a curve of marginal stability is crossed. In section five, we consider the case of an SU{2) 
gauge group in somewhat more detail. 

2. N = 2 supersymmetric gauge theories 

The field content of our model consists of an = 2 gluon vector multiplet and an 
N = 2 quark hyper multiplet. In terms of = 1 superfields, the vector multiplet may 
be decomposed as an = 1 vector multiplet V and an = 1 chiral multiplet $, both 
transforming in the adjoint of the gauge group G. The hypermultiplet consists of two 
A^ = 1 chiral multiplets Q and Q transforming in complex conjugate (in general reducible) 
representations R and R under G. To have an asymptotically free or scale invariant 
model, the index of the representation R must be less than or equal to the index of the 
adjoint representation of G. The action of the model is completely determined by A^ = 2 
supersymmetry to be of the form 

(2,1) 

where m gives the bare masses of the irreducible components of the hypermultiplet, and 
the gauge coupling constant and the theta-parameter are combined as r = ^ + Zp^. (We 
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use the conventions of for = 1 superspace. Note that our normahzation of m differs 
from that of by a factor of V2.) The classical action ( |2.1[ ) with vanishing bare masses 
is invariant under an U{l)'jz symmetry such that the 7?--characters of Wot, Q and Q 
are 1, 1/2, and respectively. The masses m and the dynamically generated scale A 
in an asymptotically free theory break this symmetry, but we may restore it by assigning 
7^-character 1 to m and A. 

The theory has a moduli space of inequivalent vacuum states. We will consider the 
Coulomb branch of this moduli space in which only the lowest component (p of the chiral 
superfield $ has a non-vanishing vacuum expectation value, subject to the constraint 
[(/), (j)^ = 0. The theory will be weakly coupled when < cj) > is large compared to the 
scale A. At a generic point on this branch, the gauge symmetry is spontaneously broken 
to its maximal abelian subgroup U{lY, where r = rank(G), and we have a vector q of 
corresponding electric charges. As usual in spontaneously broken gauge theories where 
the exact symmetry group contains abelian factors, there will also be a vector of magnetic 
charges g in the theory. In general, these charges are not well-defined over the moduli 
space of vacua; if we assemble q and g to an electric-magnetic charge vector Q = {q,g), 
this ambiguity amounts to the action of an element of Sp{2r, Z) on Q as we encircle a 
singularity 0. Finally, there is a set of conserved quark number charges S, one for each 
irreducible component of the representation R. These charges can pick up contributions 
proportional to q and g as we encircle a singularity of the moduli space . 

3. The BPS spectrum 

The central charge Z for a state is a linear combination of the electric, magnetic and 
quark number charges: 

Z = a-q + aD-g + rn-S, (3.1) 

where classically a is given by the eigenvalues of (j) and ao = Ta, but these coefficients 
get modified at the quantum level. The formula (|3.1j ) is derived by explicitly constructing 
the supercharges and calculating their Poisson bracket 0. (The contributions from 
electric and quark number charges are most easily checked by calculating the anticommu- 
tator of two supersymmetry transformations on the component fields of the vector and 
hypermultiplets. ) 

We will now determine the spectrum of BPS states in the theory, at least at weak 
coupling (i.e. for <(p> large compared to A in asymptotically free theories or Imr large) 
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where semiclassical methods can be trusted. We start with the elementary field quanta. 
The fields of the N = 2 vector multiplet fill out a BPS representation of maximum spin 
1. Indeed, the bound ( |1.2|) is saturated by the usual Higgs formula for the mass of gauge 



bosons in a situation with spontaneous symmetry breaking [|T^ . The N = 2 hypermultiplet 
must also give rise to BPS representations, since it only contains fields of maximum spin 
1/2. The central charges of the different particles are given by the eigenvalues of the matrix 



(f) + and the masses, which can be read off from the action (|2.1|) , again saturate ( p..2| ). 

The theory also contains magnetically charged solitonic excitations with non- vanishing 
expectation values of (p and the gauge field A^. The bound ( |1.2[ ) is saturated by the 



Bogomolny-Prasad-Sommerfield limit [11 1 ||T3] of monopoles and dyons |T3[ - |Tq]. Such a 



configuration breaks half of the supersymmetries Acting with the broken supersym- 
metry generators on a state of spin 0, we can create two helicity states transforming as spin 
1/2 and one more spin state, i.e. a BPS multiplet of maximum spin 1/2. We should also 
include the CPT conjugate states of opposite electric and magnetic charges for a total of 
eight helicity states. After an electric-magnetic duality transformation, these states could 
be described by a hypermultiplet coupled to the dual photons in a t/ (1)'' gauge theory 0. 

Another type of BPS states can be thought of as bound states in a quark-soliton 
system. The fermionic fields in the (Q, Q) hypermultiplet make up a Dirac spinor t/; and 
a conjugate spinor i/j = transforming in the R and R representations respectively. 

The equation of motion for these fields derived from (|2.1[) in a fixed and (p background 
reads 

{i'y^D^ - Re{(j) + m) + 7^Im((^ + m)) ip = 0, (3.2) 

where we have decomposed the matrix (p + m in its Hermitian and anti-Hermitian parts. 
A convenient basis for the Dirac matrices is 



(3.3) 

where the cr* are the Pauli matrices. We decompose the Dirac spinor accordingly as 
ip = [ I _ . Choosing the gauge Aq = 0, replacing —ido by the energy eigenvalue E and 



multiplying equation (|3.2|) by 7° from the left, we get 



1^ Lt -Im{(t> + m) - E J J ~ ' 



where the operator L is given by 

L = za^A -^Re((^ + m). (3.5) 

Each normahzable solution to ( |3.4|) doubles the number of BPS solitons, since a fermionic 
state can be either empty or occupied. 



4. The discontinuity 

It might seem difficult to find solutions to (|3.4D . However, we know that the spectrum 



of BPS states cannot change unless we cross a curve of marginal stability. To determine 
the number of bound states that appear or disappear across such a curve, we choose a 
specific point on the curve where (/> is Hermitian up to a phase, which we can rotate away 
by a U{1)ti transformation so that lm.cj) = 0. If we now have a solution to Lil)~ = 0, 
a solution of ( p.4|) would be to take tj}^ = and E = — Imm. Similarly, a solution to 
L^tp~^ = gives a solution to ( |3.4|) with tp~ = and E = Imm. Note that since the central 
charges of the quark components are given by the eigenvalues of (p + m and the central 
charge Z of a magnetically charged soliton configuration in the weak coupling limit is a 
large imaginary number, the mass difference between a quark-soliton bound state and the 
soliton state should indeed equal |Imm|. 

We have thus found a space of solutions to ( |3.4| ) of dimension dimKer(L) + 
dimKer(L^). In general, this quantity is difficult to determine, but a lower bound is 
given by the absolute value of index(L) = dimKer(L) — dimKer(L^). Now we can use the 
three-dimensional version of the Callias index theorem [O : 



index(L) = — [ Tr (UdU A dU) , (4.1) 

loTTZ J 



where U = — |Re((^-|-m) |~^Re(^-|-m) and the integral is taken over a two-sphere at spatial 
infinity. The matrix U is only well-defined if the eigenvalues of Ke{(j)+'m) are non-zero; this 
is in fact a necessary condition for L to be Fredholm so that index(L) is well-defined. As 
long as we stay away from zero eigenvalues of Re(^ -|- m), index(L) depends continuously 
on (p + m and is therefore a constant. 

However, the index can change when an eigenvalue changes sign. Indeed, the central 
charges of the quark components are given by the eigenvalues of (p + m. With Im^ = 0, 
there will be a purely imaginary central charge whenever Re(^ + m) has a zero eigenvalue. 
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This aligns with the central charge of a magnetically charged soliton in the weak coupling 
limit, and a discontinuity in the BPS spectrum is therefore possible. For Rem large 
compared to the scale A, the discontinuity will take place in the weak coupling regime, 
thus justifying our reliance on semiclassical methods. We can decompose Re((^ + m) as 

Re{(l) + m) =J2^aXax\. (4-2) 

A 

where the Xa are the real eigenvalues of Ke{(f) + m) and the xa the corresponding eigenvec- 
tors normalized up to an arbitrary phase by XaXb ~ ^ab and fulfilling the completeness 
relation '^^XaXa ~ '^^^ matrix U in is then given by 

t/ = 5^sign(A^)XAXA- (4-3) 

A 

We consider a situation where one of the eigenvalues, Aq say, changes sign. A short 
calculation gives 

lim Tr {UdU A dU) ~ lim Tr {UdU A dU) = 4dxl A dxo- (4-4) 

Aq— »0+ Aq^O — 

An important point is now that xo is in general not a globally defined function over S"^ but 
a section of the line-bundle V of eigenvectors of (p with eigenvalue —Re m. The quantity 
on the right-hand side of ( |4.4| ) is globally well-defined, though, and we can calculate the 
jump in the index as 

Aindex(L) = --1^ [ dxl A dxo = ciiV)[S% (4.5) 
Zm J 

This is of course an integer as it should be, since the first Chern class of V is an element of 
H'^{S'^, Z). In physical terms, this amounts to the Dirac quantization condition: The wave 
function of a particle with electric-magnetic charge vector Q in the presence of charges 
Q' at the origin is a section of a line-bundle over — {0}, the first Chern number of 
which equals the integer symplectic product Q ■ Q'. The number of solitonic states thus 
changes by a factor of 2^^'^ as a curve of marginal stability is crossed. This number is 
invariant under the action of Sp{2r, Z) on Q and Q' and thus well-defined on the moduli 
space of vacua. One should note that although we have considered fundamental quarks 
in a soliton background in the weak coupling regime, the result has a wider validity. For 
example, as the mass parameters of the model are turned off, the elementary quarks may 
be continuously changed to magnetically charged solitons, and curves of marginal stability 
may move from the weak coupling to the strong coupling regime 0. The discontinuities 
of the BPS spectrum are unchanged, though. It would be interesting to check the validity 
of this result in cases which can not be treated at weak coupling, and also for dyon-dyon 
bound states. 
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5. The SU{2) case 

The case of an SU (2) gauge group is of particular importance, since embedded SU{2) 
monopoles and dyons seem to be the fundamental one-particle solitons also for larger gauge 
groups [1T8[| . 

A gauge invariant parametrization of the moduli space of vacua is given by the vacuum 
expectation value of -u = |Tr((^^). By using the ^''(1)7^ symmetry, we can take u to be real 
and positive. The monopole configuration is of the form 



(j) = X ■ a (pir) Ao = A, = eijk(J^ A{r), (5.1) 

where r = [x ■ x)^/^, = r~^a;* and (f){r) is real and positive. We have (j){r) — > ^Ju 
and A{r) —>■ as r —>■ cxd. This background is invariant under rotations generated by 
Ji = Li + Si + Ti, where L^, Si and Ti are the generators of orbital angular momentum, 
spin and isospin respectively. 

With an SU (2) gauge group, asymptotic freedom or scale invariance permits us to 
consider quarks in (up to four copies of) the doublet representation or (a single copy of) 
the triplet representation. The Callias formula (|4.1|) then predicts one or two normalizable 
modes respectively, which have been explicitly constructed in the case of zero bare quark 
r9| . The two modes for a triplet quark transform as spin 1/2 under J^, leading to 



mass 



solitons in BPS representations of maximum spin 1. Doublets and triplets of SU{2) are 
the only representations that arise also when we embed a monopole in an asymptotically 
free or scale invariant theory with a larger gauge group. Triplets arise when the quarks 
are in the adjoint representation of the gauge group (i.e. in = 4 theories), but also in 
SU{n) theories with quarks in the symmetric product of two fundamental representations. 

Here we will consider the case of isodoublet quarks with non-zero bare mass in some- 
what more detail. The lowest mode should share the symmetry of the monopole back- 
ground under rotations generated by J^. This leads to the Ansatz 

(V'-)a'^ = S^^^:{r) + {x ■ a)^^i;-{r), (5.2) 

where a = 1, 2 and k = 1, 2 are spinor and SU (2) doublet indices respectively. With L as 
in ( p.5|) , the equation L'i/;~ = then amounts to 



\ ^ f -2A - (P -Rem \ / ip- 

drW-)~\-Rem -2r-^ + 2A - (P J 



(5.3) 



This system of coupled linear differential equations has two linearly independent solutions. 
However, to have a regular solution at the origin, we must impose the boundary condition 
V'^(O) = 0- For large r, this solution will generically have components along both eigen- 
vectors of the matrix in ( |5.3|) , so to have a normalizable solution we must require that 
the eigenvalues Ai = —\fu — Rem and A2 = —\fu + Rem in the r 00 limit are both 
negative. This will be the case unless 

nrn 

Re4=>l, (5.4) 



where we have stated the condition in a t/(l)7^ invariant way. For fixed m, this inequality 
determines a heart-shaped region in the finite u-plane where the quark-soliton bound state 
is absent. 

Finally, we briefly comment on the quantum numbers of the soliton states. In the 
region of moduli space where we have two such states, their quark number charges S differ 
by the quark number charge of the elementary charge, i.e. by one unit. Furthermore, for 
m = 0, the theory is invariant under quark number conjugation, which fixes the quark 
numbers of the soliton states to S = —1/2 and S = 1/2 ||19|| . For non-zero m, the quark 
number of the lowest soliton state has been calculated in [2^, and ^ as m ^ 00. At 



weak coupling, the decrease of S with increasing m is rapid. 
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